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Relativistic Mean Field Theory in the rotating frame is 
used to describe superdeformed nuclei. Nuclear currents and 
the resulting spatial components of the vector meson fields 
are fully taken into account. Identical bands in neighboring 
Rare Earth nuclei are investigated and excellent agreement 
with recent experimental data is observed. 

PACS numbers : 27.70. +q, 21.10.Hw, 21.6Q.-n, 21.60.Ev, 
21.60.Jz 

Since the experimental discovery of superdeformed 
bands in rapidly rotating nuclei many unexpected fea- 
tures of these highly excited configurations have been ob- 
served (for a recent review see Ref. [||). One of the most 
striking properties is the existence of so-called identical 
bands or twin bands, i.e. nearly identical transition en- 
ergies Ej of the emitted 7-radiation in bands belonging 
to neighboring nuclei with different mass numbers. In a 
considerable number of nuclei in the Dy-region as well as 
in the Hg region one has found differences in E.^ of only 1- 
3 keV, i.e. there exists sequences of bands in neighboring 



nuclei, which are virtually identical, AEj/Ej 
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Since these transition energies are directly related to the 
corresponding dynamical moments of inertia obeying on 
the average a simple A^^^ dependence, one would have 
expected changes of one order of magnitude larger. 

Several groups have tried to understand this phe- 
nomenon by means of conventional investigations within 
the framework of the semi-phenomenological Strutinski 
method in connection with a rotating Nilsson or Saxon- 
Woods potential ^,||. It has been pointed out, that 
the different single particle orbits can give rather dif- 
ferent contributions to the moment of inertia. This is 
most clearly seen in the simple oscillator model, where 
orbits with oscillator quanta along the rotational axis 
have vanishing angular momentum operator matrix el- 
ements. Other groups realized that a new coupling 
scheme in nuclei exists, the so-called pseudo-spin scheme, 
which to a very large extent decouples the pseudo-orbital 
motion from the pseudo-spin degrees of freedom and fa- 
vors the strong coupling limit. In all these investigations, 
however, polarization effects, which are expected to pro- 
duce much larger changes of the moments of inertia than 
those observed in identical bands, are either neglected 
completely Q] or taken into account only partially by 
minimizing the rotating energy surface with respect to a 
few deformation parameters. 

We therefore feel that it is very important to carry 



out fully self-consistent microscopic calculations, where 
all the degrees of freedom are taken into account. Such 
calculations are not simple, but they are nowadays fea- 
sible. A first investigation of this type using the density 
dependent Skyrme III force has very recently been car- 
ried out in the Hg region It has been found that 
full self-consistency has indeed a considerable influence 
on the details of the moment of inertia in neighboring 
nuclei. There are, however, many important questions 
still open. In these nuclei identical bands occur in a 
spin region, where pairing plays an important role. Since 
particle-particle correlations are not very well described 
within the Skyrme scheme, it is not at all clear, to what 
extent the present deviations of the theoretical from the 
experimental results can be understood by such deficien- 
cies. 

The present investigation is therefore devoted to the 
Dy-region, where superdeformed bands are observed up 
to very high angular momenta. At these very large rota- 
tional frequencies the Coriolis-anti-pairing effect reduces 
these correlations considerably, such that they have little 
influence on the moment of inertia. We use a relativis- 
tic field theory which includes a-, to- and p-mesons as 
well as the electromagnetic field. In addition we consider 
a nonlinear self-coupling of the cr-field. We use the pa- 
rameter set NLl, which has been adjusted |6) to nuclear 
matter and a few spherical nuclei. This parameter set 
has turned out to be very successful for the description 
of many groundstate properties over the entire periodic 
table . In particular one has found excellent agreement 
with ground state deformations in open shell nuclei. The 
rotation is treated within the cranking approach, in ac- 
cordance with the concept of a mean field description. 
This leads us to a relativistic self-consistent cranking the- 
ory (RSCC) as developed in Ref. [|. 

In the rotating frame time reversal invariance is bro- 
ken. This lead to nucleonic currents in the interior of 
the nucleus, which form the source of magnetic poten- 
tials in the Dirac equation (nuclear magnetism). In this 
way the charge current jc is the source of the normal 
magnetic potential A, the isoscalar baryon current js is 
the source of the spatial components uj of the w-mesons 
and the isovector baryon current ja is the source of the 
spatial component ps of the p-mesons. In contrast to the 
maxwellian magnetic field A having a small electromag- 
netic coupling, the large coupling constants of the strong 
interaction causes the fields <jJ and p to be important in 
all cases, where they are not forbidden by symmetries. 



1 



such as time reversal. They have a strong influence on 
the magnetic moments in odd mass nuclei, where time 
reversal is broken by the odd particle, as well as the mo- 
ment of inertia in rotating nuclei, where time reversal 
is broken by the Coriolis field. In an early investiga- 
tion of rapidly rotating superdeformed nuclei within the 
framework of cranked relativistic mean field theory @ 
these components were not taken into account for rea- 
sons of numerical simplicity. Strong deviations from the 
experimentally observed moments of inertia were found. 
Only the size of the quadruple deformation was repro- 
duced properly. Semiclassical corrections turned out to 
be large, but could not reproduce the proper experimen- 
tal values of the moment of inertia. 

In this investigation nuclear magnetism is taken fully 
into account in a self-consistent way. Starting from the 
Lagrangian 

C = -il:(^- g^tp + gppf ~ ie(l + Tg)^ - g„a - Mj^ ^ 
+ \d^(Td^'cj - U{a) - if^^.r!^'^ + ^m^c^^c^^ (1) 

where Mn is the bare nucleon mass and '0 is its Dirac 
spinor. We have in addition the scalar meson (ct), 
isoscalar vector mesons [lo^), isovector vector mesons 
{f?^) and the photons A'', with the masses mo-, rn^ and 
rUp and the coupling constants g^, gui, gp- For simplicity 
in the following equations we neglect the p-meson and 
the photon. In the calculations these contributions are, 
however, taken into account however. The field tensors 
for the vector mesons are given as 



(2) 



For a realistic description of nuclear properties a nonlin- 
ear self-coupling for the scalar mesons has turned out to 
be crucial |13|: 
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(3) 



Starting from this Lagrangian and transforming to a 
frame rotating with a uniform velocity around the x- 
axis perpendicular to the symmetry axis of the deformed 
nucleus in its ground state, we obtain the classical equa- 
tions of motion 

{a(p -I- g^(jj) + g^ujQ + (3{M + g^a) - n^J^jipa = ^a'lpa 

(4) 

for the nucleon spinors and 

{-A + {n,Lxf}a + U'ia) 
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(5a) 
(5b) 
(5c) 



where = + Sx and the spin operator S'^, is a 4 x 4- 
matrix for the spinor fields with spin i and a 3 x 3-matrix 
for vector fields with spin 1. For details see Ref. 

These equations are solved self-consistently by expand- 
ing the Dirac spinors as well as the meson fields in terms 
of eigenfunctions of a deformed oscillator, as discussed in 
details in Refs. Up to Np — 13 major oscillator 

shells where taken into account for the large components 
of the Fermions fields and up to Nb = 13 shells for the 
meson fields. Because of the large number of configura- 
tions in this spaces high-lying orbitals with a deformed 
oscillator energy larger than 10.3 x huj'^ for the Fermions 
and larger than 10.5 x huj'^ for the Bosons have been ne- 
glected. In general one has to allow for complex expan- 
sion coefficients in this basis. However, assuming mirror 
symmetry at the three planes {x,y), (a;, z), and {y,z) 
for the densities and rotational symmetry for the cur- 
rents it is possible to restrict oneself to real coefficients. 
Since, in the rotating frame, their is no a priori reason for 
these symmetries, a complex code was used to show that, 
even with initial conditions strongly violating this sym- 
metries, we find after many iterations final self-consistent 
solutions which obey these symmetries. All the following 
calculations have therefore been carried out with the real 
code only. 

In Fig. 1 we show the static and the dynamic mo- 
ment of inertia for the lowest superdeformed band in 
the nucleus ^^^Dy as a function of the angular momen- 
tum. It is clearly seen that a calculation without nu- 
clear magnetism, i.e. without the spatial contributions of 
the vector meson-fields, which is in good agreement with 
the experimental quadrupole moments (see Ref. [|o|) 
produces much too small moments of inertia. A semi- 
classical correction where these contributions, derived in 
Thomas Fermi approximation using a rigid rotor current, 
are taken into account in first order perturbation theory, 
overemphasized the moments of inertia by roughly 10 
%. Only if one takes these contributions into account 
in a fully self-consistent way, is perfect agreement with 
experimental data achieved. In the region of small an- 
gular momenta one still observes very small deviations, 
which could possibly be understood as the influence of 
remaining pairing correlations in this region of interme- 
diate spins. 

We also find that nuclear magnetism practically has 
no infiuence on the shape of the nucleus. The mass 
quadrupole moments decrease in the spin range from 20 
to 60 h only very little, running from 4350 to 4260 fm^ 
and the corresponding hexadecupole moments change for 
the same region from 20600 to 19600 fm"*. The changes 
induced by nuclear magnetism are of the order of a few 
per mille. The average charge quadrupole moment if 
found to be 18.6 eb, which is in good agreement with 
a value of 18 eb obtained in a non-relativistic calcula- 
tion jll], and the experimental value of 19 eb |jl^ From 
the quadrupole moments we can derive the Hill- Wheeler 
parameters (3 = 0.72 and 7 = 0.7° for the quadrupole de- 
formations, which corresponds closely to a nearly prolate 
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deformed nucleus with an axis ratio of 1:1.9, close to the 
standard value 1:2 of the harmonic oscillator model. 

Let us now investigate the problem of identical bands. 
For this purpose we calculate, in a self-consistent way, 
bands in the neighboring nucleus ^^^Tb by removing one 
proton from the ^^^Dy core. In Fig. 2 we show the single 
particle spectrum for protons in the rotating potential 
formed by the ^^^Dy core. The large gap at Z = 66 
is clearly recognized. Taking particles out of the orbits 
directly below this gap, we can produce different bands 
in the nucleus ^^^Tb. They have the quantum numbers 
(PS) of parity (P) and signature (S), namely (H — ) for the 
dashed line, (++) for the full line, ( — h) for the dotted 
line and {—+) for the dashed dotted line. 

The proton hole induces a polarization of the ^^^Dy- 
core, which has two effects: it leads to changes of defor- 
mation and in addition to changes in the current distribu- 
tion. In Table | we show the values obtained after solving 
in a fully self-consistent fashion the relativistic mean field 
equations for the odd system in the four lowest configu- 
rations. We show the values for several observables for 
the lowest superdeformed band in ^^■^Dy. The relative 
changes with respect to this reference band in the four 
bands of ^^^Tb are given in per mille. According to the 
simple A^^^-Tule we expect changes in the moment of in- 
ertia by « 11 per mille. The calculated values for the 
moments of inertia J'^ and for the band with the 
quantum numbers ( — h), which we we shall in the fol- 
lowing call the identical band, are, however, at least an 
order of magnitude smaller. This is by no means triv- 
ial, because we find considerably larger changes in the 
quadrupole moments and in the rigid body moments of 
inertia. In fact in most of the other bands the changes 
are also much larger. 

In order to have a direct comparison with the ex- 
periment we show in Fig. 3 the differences AE-y — 

(Tb) — E^ (Dy) between the transition energies in sev- 
eral bands in the nucleus ^^^Tb and in the lowest superde- 
formed band in ^'''^Dy. The agreement with the experi- 
mental value is excellent for the band with the quantum 
numbers ( — h), where the energy differences are of order 
of 1 keV. As we see in Fig. 2, this band correspond to a 
hole in the orbit with the approximate Nilsson quantum 
numbers [301]^ — . This orbit has a very small number of 
oscillator quanta along the z-axis (the symmetry axis), 
which yields nearly vanishing contributions to the mo- 
ment of inertia. We are therefore in agreement with the 
qualitative argument put forward in Ref. fl^ 

This is, however, not the full story. In order to in- 
vestigate the very good quantitative agreement, we have 
carried out two additional calculations in Fig. 4 for the 
identical band band with the quantum numbers ( — h): 
First we neglected the polarization induced by the pro- 
ton hole, i.e we calculated the energy differences for wave 
functions for the nucleus ^^^Tb obtained from the ^^^Dy 
core by just removing one proton, without requiring self- 
consistency for the odd mass configuration. In this case 
we find the dashed dotted line in Fig. 4, which is in sharp 



disagreement with the experimental data. Next we took 
into account the polarization, but we neglected nuclear 
magnetism, i.e. the contributions of the spatial compo- 
nents of the vector meson fields and find the dashed line 
in Fig. 4, which is also in disagreement with experiment. 

We therefore conclude, that a very delicate cancella- 
tion process occurs in identical bands in superdeformed 
nuclei. Polarization of the quadrupole moments and of 
the density alone would induce changes of the order of 
5-10 per mille. Neglecting nuclear magnetism would 
also lead to changes of this order of magnitude. Obvi- 
ously both act in opposite direction, such that the finial 
differences are only in the order of 1 per mille. So far 
the precise mechanism for this cancellation is not fully 
understood. It requires definitely much more systematic 
investigations. Nonetheless it seems to us a very satisfy- 
ing and surprising result, that without any free parame- 
ter, and simmply using the set NLl adjusted to nuclear 
matter and a few spherical nuclei, long before identical 
bands had been identified, we can obtain this degree of 
accuracy in the relatively simple minded relativistic mean 
field approach. We have to emphasize, however, that full 
self-consistency as well as the inclusion of the nuclear 
currents are very important in this context. 
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Figure Captions 



Table Caption 



Fig. 1 (a) Static (J^^)) and (b) dynamic (J^^)) 

moment 

of inertia for the lowest superdeformed band in the 
nucleus ^^^Dy. The dashed line corresponds to the 
calculation without the spatial contributions of the 
vector mesons. In the dotted line such contribu- 
tions are taken into account in a semiclassical way 
and the full line represents the self-consistent solu- 
tion including these contributions fully. 

Fig. 2 Single proton spectra in the self-consistent rotat- 
ing potential of the superdeformed band in ^^^Dy as 
a function of the cranking frequency Clx- For small 
values of fl^, where the two signatures are nearly 
degenerate, we indicate the approximate Nilsson 
quantum numbers. 

Fig. 3 Differences AE^ in the transitional energies — 
£'(/)—£'(/ — 1) between the lowest bands for each 
pair of quantum numbers (P,S) in the nucleus ^^^Tb 
and the superdeformed band in ^^^Dy are compared 
with experimental values for the excited superde- 
formed band with negative parity in ^^^Tb. 

Fig. 4 Differences AE^ for the identical band with the 

quantum numbers (— +). The fully self-consistent 
solution (full line) and solutions neglecting nuclear 
magnetism (dashed line) or polarization induced by 
the proton hole (dashed dotted line) are compared 
with the experiment (empty diamonds). 



TABLE I. Binding energy E, mass quadrupole moment 
Qo, dynamic (i7^), static {J^^) and rigid body {Jrig) moment 
of inertia at the angular momentum / = 507i for the superde- 
formed band in ^^'^Dy and relative changes of these values 
(given in %) several bands in the neighboring nucleus ^^^Tb 
at the same angular momentum. 



Band E (MeV) Qo (fm^) J^^' J^^' J^, 



152 Dy 


-1228.358 


4287.7 


82.544 86.41 


93.35 


i5iTb*(-|-,-|-) 


0.53 


-2.91 


-1.554 -1.25 


-1.69 


i^^Tb (+,-) 


0.51 


-3.17 


0.654 0.35 


-1.80 


151-pb*(_,+) 


0.56 


1.30 


-0.001 0.10 


-0.45 


i5iTb*(-,-) 


0.54 


1.19 


0.145 1.48 


0.39 
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